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Abstract 

A new scaling variable is introduced in terms of which nuclear shadowing in deep- 
inelastic scattering is universal, i.e. independent of A, and x. This variable can be 
interpreted as a measure of the number of gluons probed by the hadronic fluctuations 
of a virtual photon during their lifetime. The shadowing correction grows at small x 
substantially less steeply than is suggested by the eikonal approximation. This results 
from the fact that shadowing is dominated by soft rather than hard interactions. 
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1. Introduction 

A large amount of new high precision data on nuclear shadowing in deep-inelastic scattering 
is now available. The comparison with theory is not easy since the data are taken for a 
variety of nuclei at different values of and x. One of the goals of this paper is to find a 
scaling variable, dependent on Q^, x and A, which makes shadowing an universal function 
of this variable. Our guess of the correlations between Q^, x and A entering the scaling 
variable is inspired by the prejudice about the underlying QCD dynamics of shadowing. 

Deep-inelastic scattering is usually interpreted in two alternative frames of reference, in 
the infinite momentum frame of the proton in terms of its structure function, or, in the 
proton rest frame in terms of hadronic fluctuations of the photon. In the first case, nuclear 
shadowing looks like a result of overlap in the longitudinal direction of the parton clouds 
originated from different bound nucleons. The latter approach is more in line with the 
familiar vector dominance model It seems to be advantageous treating shadowing at 
small X as the consequence of the hadronic fluctuations of the photon. On the technical 
level our approach is similar to that of 0. 

Prior to the detailed discussion, we introduce the scaling variable motivated by the 
Glauber model as an average number of bound nucleons taking part in interacting with 
the hadronic fluctuation of the photon, 

n{x,Q\A) = \^{T^)Fl{q,) (1) 

In this expression is the total cross section of interaction of the hadronic fluctuation 
h of the photon with a nucleon. We choose the basis of states h to be eigenstates of the 
interaction, rather than of the mass matrix. This eliminates the off diagonal amplitudes in 
the double dispersion relation |^ for deep-inelastic scattering. The concrete choice of the 
eigenstate basis is done below. The averaging in eq. (|I]) is weighted with the probability 
to find the fiuctuation h in the photon. (T4) is the mean nuclear thickness function of the 
nucleus. The nuclear formfactor -Fa(q'l) depends on the longitudinal momentum transfer in 
the diffractive dissociation 'jNQhN. 

In what follows, we try to justify our choice of the scaling variable and provide model- 
dependent estimates of n{x, Q'^,A). 
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2. A model for n(x,Q'^,A) 
At high energies the hfetime of the hadronic fluctuations of the photon may substantially 
exceed the nuclear radius, 2z//Q^ ^ Ra- In this case, the nuclear photoabsorption cross 
section can be represented in the same form as in the hadron-nucleus interaction ^, ^, 



Here T{b) ^ dzpA^b, z) is the nuclear thickness function, where z) is the nuclear 
density, which depends on the impact parameter h and the longitudinal coordinate z. The 
dipole cross section cr(p, x) of the interaction with a nucleon of a qq pair depends on its 
transverse separation p and the energy related to a; = jlvnjqv . The averaging over p and 
a, the fraction of the photon light-cone momentum carried by the quark, is weighted with 
the photon wave function squared. One should be cautious with such a definition of the 
averaging [§] in the case of the photon, because its wave function, being different from the 
hadronic one, is not normalized to one. At small p it has a form |\E'^, (p, = 
6aem/(2vr)^Sjej[l — 2a(l — a)] i^K\{ep\ where = a{\ — a)Q'^ \rr?^. The mean transverse 
separation of the qq fiuctuation is given by (p^) oc and is of the order l/Q"^. However, 
at the endpoints of the kinematical region a or 1 — a ~ mg/Q^, the qq fiuctuations acquire a 
large transverse size, p^ ~ V""^? 0) S H- For light quarks such a big size may substantially 
exceed the confinement radius, and one should put a cut off on the integration over a. This 
is equivalent to a replacement of the quark mass by the cut off A. 

Expanding eq. @, we represent the nuclear photoabsorption cross section in the form 



tot {^^ Q"^) 



2 d%il- 



a{p,x)T{b) 
2A 



^L^ix, Q') = A a7o7(x, Q') [l - nix, Q^ A) +...], (3) 

where n{x,Q\A) = {T{b)){a^) /A{a) and (T(6)) = {A - 1)/A^ J d%T'^{b). 

Note that expansion (|^) looks similar to that given by the Glauber approximation 0]. 
In fact we use the eikonal Glauber formalism for projectile states with definite transverse 
dimension p, since they are the eigenstates of interaction. However, as we conclude below, 
after averaging over the photon wave function the first shadowing correction is of the order 
of rather than l/Q"^ as expected in the Glauber model. This comparison demonstrates 
that, in terms of multiple scattering theory, DIS on nuclei is dominated by Gribov's inelastic 



shadowing 0, while the Glauber eikonal contribution vanishes at high Q"^. In such a 
case , A) should be interpreted as a measure of a number of gluons probed by the qq 

fluctuation of the photon, rather than a number of nucleons. This is justified at small p since 
the photoabsorption cross section is proportional to and the gluon distribution function 



10| , PJ| . The latter is not proportional to the nucleon density because of gluon fusion [|12 



- [T^. This effect is related to the inelastic corrections corresponding to the excitation of 
heavy mass intermediate states. On the other hand, if p is not small, the qq fluctuation 
experiences additional shadowing interacting with gluons. This results in a rather small 
unitarity correction to the photoabsorption cross section on a nucleon, but in a substantial 
correction on nuclei. Taking into account the shadowing corrections coming from large size 
fluctuations, we may say that such qq fluctuations also probe the number of gluons. 

In order to evaluate the lowest order nuclear shadowing correction n{x^ Q^, A) in eq. (^), 
note that the denominator in eq. is directly related to the proton structure function, 
{a{x,p)) = 'in'^aem F2{x,Q^)/Q'^. We performed a fit to available data on F^a;, Q^) from 
NMC [0, HI [T^, ig and ZEUS 0, ^ experiments with x < 0.05 and > O.b/GeV^. We 
used a simple parameterization motivated by the double-leading-log approximation (DLLA) 
for QCD evolution equations [|1^] (see also review 0]), F^x, Q^) = f{Q'^)[aexp{2\^)/L + 
b^/x\, where the first term corresponds to the sea- quark (Pomeron), while the second term 
originates from the valence - quark (Reggeons) contributions. L = {Att/Pq) \n{c/as) \n{d/x), 
where = 47r//3o ln(Q^/AQ(^^) and Aqcd = 0.2 GeV, Po = 9 for three active flavors. The 
factor f{Q'^) = Q'^/{e + Q"^) guarantees that the structure function vanishes in the limit of 
real photoabsorption. This parameterization fits the data very well {x^/d.f. = 0.8) with 
parameters a = 0.036 ± 0.005 6 = 0.4 ± 0.08, c = 0.59 ± 0.02, d = 0.31 ± 0.07, e = 
0.12 ±0.08 GeV^. 

Although DLLA does not provide a Regge form of the structure functions, an effective 
Regge parameterization may be a good approximation, F2{x,Q^) oc exp[Ae//((5^)^]. Here 
^ = ln(l/x) and Aeff{Q'^) = dln[F2{x,Q'^)]/d^ correspond to the effective Pomeron inter- 
cept tteff = 1 + Ag//. Of course, Ae// can be treated as x- independent only in a restricted 
interval of x. The values of Af.ff{Q'^,x) corresponding to the results of our fit are depicted 
in fig.l as function of versus x. We see that Ae//((5^,x) is almost x-independent which 
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justifies the Regge parameterization as a good approximation in the range of x and Q'^ under 
consideration. 



0.4 




Figure 1: The effective Pomeron intercept as function of versus 
X = 10"^'^'"^'^. The proton structure function (x, Q^) is fitted to avail- 
able data as is explained in the text. 



Remarkably, the values of Ag// in fig. 1 substantially exceed what is known from the 
energy dependence of the total cross sections of proton-proton interaction |2^, |2^, Ag// ~ 



0.07 — 0.08. Our results show that this distinction remains substantial down to quite low 
~ 2 GeV^. Note that the observed Q^-dependence of Ae// contradicts the Pomeron 



factorization. This is not surprising since perturbative QCD calculations ||2^, ^ show that 
the Pomeron is a more complicated singularity, a cut or a sequence of poles. 

The rising Q^-dependence of Ae//(Q^) means that the x-dependence of the cross section 
a{p,x) is steeper at smaller p. Indeed, the larger Q"^ is, the smaller is (p^) oc ln(Q^/A^)/Q^. 

Despite the smallness of the mean transverse size of the photon fluctuations partici- 
pating in DIS, of the order of ~ the shadowing terms in the expansion eq. (|^) are 
dominated by large transverse separations in qq fluctuations. This can be argued using 
the relation, (p^)/(p^) = 2.4/A^ln((5^/A^), which is obtained by the same averaging pro- 
cedure as is defined in eq. (^. This relation is the manifestation of a salient feature of 
the hadronic fluctuations of the photon, namely, a huge dispersion of the transverse size 



distribution, (p^) ^ (p^)^- Although we used the perturbative photon wave function, this 
conclusion has a rather general character. It is a result of the interplay of perturbative 
and nonperturbative contributions in the deep-inelastic cross section. The former results 
from the small-size fluctuations corresponding to "symmetric pairs", a ~ 1 — a. They are 
presented with large probability in the photon wave function, but have a small, ~ l/<5^ 
interaction cross section. On the contrary, the highly asymmetric fluctuations with a or 
(1 — a) ~ A^/Q^, have a very small, ~ weight in the photon wave function, but a large 
interaction cross section, typical for hadrons. In the case of double scattering in the nucleus 
the perturbative contribution turns out to be very much suppressed by the factor of l/Q"^, 
while the nonperturbative part is suppressed only once by the weight factor of l/Q^. Thus, 
the soft interaction dominates nuclear shadowing. 

Actually, just this effect is responsible for the scaling behavior of nuclear shadowing p 
and of unitarity corrections to the photoabsorption cross section on a nucleonQ. 

Once the interaction responsible for shadowing is essentially soft, its x-dependence is 
governed by the soft Aejj(A^) ~ 0.1, rather than the hard one. This is confirmed by 
the recent study of diffractive dissociation by the HI collaboration P^, which claimed 
Ap = 0.1 ±0.03 ±0.04. 

To proceed further with the calculation of n{x,Q'^,A), note that in eqs. (0), ® we 
temporarily used an assumption that the photon energy u = Q'^/2m]\fX in the nuclear rest 
frame is sufficiently high to make the lifetime of the photon fluctuation long compared with 
the nuclear size, so that it propagates through the whole nucleus with a frozen intrinsic 
separation p. However, most of the data available are in the transition region of x, where 
the lifetime, usually called coherence time, is comparable with the nuclear radius. The 
finite coherence time can be taken into account by introducing a phase shift between qq 
wave packets produced at different longitudinal coordinates, in the same way as for inelastic 
corrections p8| , or in the vector dominance model [jl| . The mean nuclear thickness function 
of eq. (H) should be replaced by an effective one. 



■fThis conclusion is in variance with the statement in p6|| that at high the unitarity corrections vanish 
as ^ l/Q^ and one sees the single Pomeron exchange. 
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' ^ (T(6)) Fl{q,) . (4) 

Here FA^qh) = 6xp(— g|_R^/6) is the nuclear longitudinal formfactor and Ra is the mean 
nuclear radius. For the sake of simplicity we use the Gaussian form for the nuclear density 
which is quite precise for the calculation of FA^qi)- Calculating (T) we use the realistic 
parameterization of nuclear density p9| . 

The decrease of the effective nuclear thickness function {T{h)) at large qi can be inter- 
preted as a result of shortness of the hadronic fluctuation path in the nucleus, if we are in 
the nuclear rest frame, or as an incomplete overlap of the gluon clouds of the nucleons which 
have the same impact parameter in the infinite momentum frame of the nucleus |l30| , ^ . 

In order to calculate the longitudinal momentum transfer in DIS, = (Q^ + M'^)/2v, 
one needs to know the effective mass of the produced qq wave packet. However, a qq state 
with definite separation p does not have a definite mass. This is a typical problem for 
those who work in the eigenstate basis of interaction. We evaluate qi ~ 2xmN assuming 
~ Q^. Thus, the parameter which controls the value of {T{b)) is x. 

Now we are in a position to estimate the nuclear shadowing correction per nucleon 
n{x,Q'^,A) in eq. @, which reads 

1 AT /1\2Ae//(A2) 

The scaling variable n{x, Q^, A) can be interpreted as a measure of the amount of those 
gluons which take part in the interaction with the qq fluctuation during its lifetime. Eq. (|^) 
is a model-dependent realization of the general expression (p. 

3. Comparison with the data 

The variable n{x, Q^, A) has been calculated from eq. (P) using the values of x, and A 
corresponding to data from the NMC experiment Q as well as the results of our fit to 
(x, Q^). The expected scaling dependence of the nuclear shadowing on n{x, Q^, A) is not 
affected by the overall normalization factor N. However, in order to have a correct slope 
of n-dependence corresponding to eq. @ we choose = 3 GeV~^. The data on the ratio 
of the nuclear and nucleonic photoabsorption cross sections Ra/n{x,Q^) is plotted against 



A - i 

~A^ 



dz pA^b, z) e 



iqz 
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n{x,Q'^,A) in Fig. 2. They demonstrate a good scaling in n{x,Q'^,A) within a few percent 
accuracy. 
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Figure 2: Data on nuclear shadowing at small x from the NMC ^ 
experiments versus the scaling variable n{x, Q^, A) as defined in eq. (^. 
The slope of the straight line corresponds to the normalization factor in 
eq. d) AT = 3 GeV-\ 



We should comment more on the procedure of the calculation of n{x, Q^, A): 
(i) Our considerations are valid only for small x, so we limit the x-region to x < 0.07. 
At larger x, the nuclear structure functions show a small enhancement of a few percent 
relative to the proton one, which results in > 1 for Q^, ^)fiO. A plausible 

assumption is that about the same antishadowing correction extends down to smaller x, 
where it is compensated by stronger shadowing effects. Such a behavior, for instance, is 



expected in the model of swelling bound nucleons |3J, |10|- The antishadowing effect may 
have some A dependence, what would cause a small, a few percent relative shift of the data 
in fig. 2 corresponding to different nuclei, but will not change the slope of n-dependence. 
Since the physics of antishadowing is beyond the scope of our present consideration, and 
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the effect is numerically very small, we do not try to incorporate with it, but just have 
renormalized the solid line Ra/nIu) = 1 — n in fig. 2 by 3% up to make the comparison 
easier. 

(ii) The data points H for < 0.5 GeV^ were excluded from the analysis because they 



are in the realm of the vector dominance model, rather than DIS. They should correspond 
to the same nuclear shadowing as is experienced by the p-meson. This is the reason for the 
saturation of nuclear shadowing at small x, claimed in ^ . 

(iii) A further important observation is that Ra/n{x,Q'^) depends to a good accuracy 
linearly upon n{x,Q'^,A) at least for n < 0.2. The higher order terms in the expansion in 
eq. @ are expected to violate the linearity in n. This implies that those terms are small. A 
model-dependent evaluation of the next shadowing correction shows that it is small indeed. 



(iv) There are two contributions to the shadowing in DIS ||T2|, |13|, one comes from the 
suppression of the gluon density as a consequence of gluon fusion ggQg, which corresponds 
to the triple Pomeron graph in the framework of standard Regge phenomenology. Another 
contribution to the shadowing comes from the Glauber-like rescattering of the qq fluctuation 
off gluons. The latter mechanism, which was mostly under consideration above, can also be 
viewed upon as a parton fusion, but as a fusion of gluons into a qq pair, ggQqq. In the Regge- 
model this process corresponds to the Pomeron- Pomeron- Reggeon graph. Both mechanisms 
lead to the same form of the variable n{x, Q^, A) in eq. (^). However, the formfactor Ffyiqi) 
has a different form for the triple-Pomeron mechanism due to the contribution of heavier 
hadronic fluctuations of the photon, FAiqi) oc Ei{—x'^m%R\/6). Here Ei is the integral 
exponential function. We checked that in the x and A domain investigated, the admixture of 
the triple-Pomeron mechanism does not affect the n{x, Q^, A)-scaling within the error bars 
of the data available. It may, however, cause a deviation from the scaling for heavy nuclei. 
We hope that forthcoming high-statistic data on heavy nuclei from the NMC Collaboration 
may help to disentangle these two mechanisms of shadowing. This is important if one wants 
to predict the unitarity corrections to the proton structure function at small x or the photon 
diffractive dissociation cross section, because the admixture of the triple-Pomeron affects 
the normalization constant in eq. (^. 

4. Summary 
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Starting from the QCD dynamics of deep-inelastic scattering at small x we have found a new 
variable n{x, Q^, A) which scales all available data on nuclear shadowing in DIS at small-x. 
This variable measures the number of gluons which a hadronic fluctuation of the virtual 
photon interacts with during its lifetime. 

An important observation is also that shadowing corrections at small x and large Q'^ 
grow less steeply than [FK^, Q'^Yf- This is because nuclear shadowing is a subject to soft 
rather than hard physics. 

The observed scahng of nuclear shadowing as function of n{x,Q'^,A) supports our as- 
sumptions on the dynamics of nuclear shadowing. 
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